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MODAL DENSITY ESTIMATES FOR SANDWICH PANELS: 


THEORY AND EXPERIMENT 1 
Larry L. Erickson 
Ames Research Center 


SUMMARY 


Formulas are presented that can be used to estimate the average modal 
densities of sandwich beams and flat or cylindrical ly curved sandwich panels. 
Numerical results, presented in terms of general parameters, indicate the rela- 
tive importance of transverse shear flexibility, orthotropic shear moduli of 
the core, face bending stiffness, rotary inertia, and panel curvature over a 
wide frequency range. Modal densities determined experimentally from resonance 
tests of flat rectangular sandwich panels having orthotropic cores are close to 
the modal densities estimated by theory except at frequencies near that of the 
fundamental mode. In this low-frequency range the theoretical estimates are 
not valid. 


INTRODUCTION 


Predicting the vibratory response of structures to random environmental 
loads that are distributed over a wide frequency range (e.g., rocket engine 
and turbulent boundary- layer noise) presents considerable difficulties. In 
the conventional approach to vibration analysis the response is expressed as a 
series, the terms of which involve the structure's free vibration mode shapes 
and resonance frequencies. Unfortunately, for most practical structures the 
mode shapes and resonance frequencies, especially those occurring at the 
higher frequencies, are difficult to obtain. As a consequence, there has been 
an effort during the past few years to develop a new approach to multimodal 
vibration problems that avoids the difficulties of expanding the response in 
terms of the mode shapes (see, e.g., refs. 1 and 2). In this approach, some- 
times referred to as "statistical energy analysis," average response levels in 
various frequency bands are estimated without knowledge of the mode shapes and 
resonance frequencies. Instead, what is required is a knowledge of the type 
(e.g., flexural) and number of structural vibration modes occurring in a given 
frequency interval. The number of modes per unit frequency is called the 
"modal density" of the structure. 

Because a structure's modal density is relatively independent of the 
boundary conditions, statistical ene rgy analysis shows promise of becoming a 

*Part of the information presented herein was published under the title 
"Modal Densities of Sandwich Panels: Theory and Experiment," in the Shock and 

Vibration Bulletin, Bull. 39, Pt. 3, U.S. Dept, of Defense, Jan. 1969, 

pp. 1-16. 



useful analysis technique for estimating average response levels of multimodal 
structural vibrations. This approach has been used to study a number of prac- 
tical vibration problems including launch vehicle response to acoustic pres- 
sure fields (refs. 3 and 4), noise vibration transmission in space vehicles 
(ref. 5) , vibration transmission into an instrument package (ref. 6) , and 
sound transmission through a partition (ref. 7). In each case, the modal den- 
sities of the structures were required in the analysis. For example, one of 
the more useful response quantities is acceleration spectral density, which is 
directly proportional to modal density. 

Equations have been derived for the modal densities of several structural 
elements and can be found in references 8 and 9. A portion of the results in 
reference 9 pertains to flat and doubly curved sandwich panels, and it is 
shown there that the transverse shear flexibility of a panel can have a signif- 
icant effect on the panel's modal density. However, the results of refer- 
ence 9 apply only to sandwich panels with isotropic cores and whose faces 
behave as membranes (i.e., the face bending stiffness is neglected). Since 
the cores of lightweight sandwich panels often have shear orthotropy (e.g., 
honeycomb) , there is a practical need for equations that account for the 
effect of core orthotropy on the modal density. In addition, it is known from 
sandwich-panel buckling theory that for small wavelengths of deformation the 
face bending stiffness can become important (ref. 10) . 

In the present investigation, modal density estimates are obtained analyt- 
ically for sandwich beams and for flat or cylindrical ly curved sandwich panels. 
The effects of the core's transverse shear flexibility, including orthotropic 
shear moduli parallel to the faces, face bending stiffness, and rotary inertia, 
are examined. Numerical results are presented in graphic form. 

In addition, experimental results obtained from resonance tests of flat 
rectangular sandwich panels having orthotropic honeycomb cores are presented. 

In these tests, up to 80 consecutive vibration modes per panel were excited 
and identified. The experimentally determined modal densities are compared 
with the modal densities predicted by theory. 


THEORY AND ASSUMPTIONS 


Panel 

The structural elements considered are a flat or cylindrical ly curved 
sandwich panel and a sandwich beam. The coordinate system and geometry of the 
curved sandwich panel are shown in figure 1. The panel is of length a and 
circumferential width b, and has a constant radius of curvature R. Each 
face sheet is isotropic, homogeneous, and of uniform thickness. The thickness. 
Young's modulus, and mass density of the top face sheet may be different from 
the thickness, modulus, and density of the bottom face sheet. The core may 
possess, orthotropic shear moduli (G Cx , G C y) and is of uniform thickness h c> 


r 


The modal density estimates obtained herein are based on three frequency 
equations that relate the panel natural frequencies (co rad/sec) to various 
physical characteristics of the panel and to the number of half-waves m and n 
that form in the x and y directions, respectively. The first of these 
equations, which is derived from the theory of reference 11 (see appendix B) , 
is 



where p 2 = [m/(a/b)] 2 + n 2 . In this equation, the physical and geometric 
properties of the panel are described by the following dimensionless 
quantities : 

1. The length-width ratio a/b 

2. The curvature parameter C, which is proportional to 1/R 

3. The rotary inertia parameter X, which is proportional to the panel's 
mass density moment of inertia 

4. The shear flexibility parameters r x and ry, which are proportional 
to l/G Cx and 1/G C ^, respectively 

The reference frequency 0 ) o is the fundamental frequency of a flat, 
simply supported, semi-infinite sandwich panel with length b (a/b = °°) as 
predicted by classical plate theory (shear and rotary inertia neglected); p is 
Poisson's ratio of the faces and is assumed to be the same for both faces. 

In equation (1), the bending stiffness of the faces about their own middle 
surface has been neglected. The second frequency equation accounts for the 
face bending stiffness but does not incorporate the effects of rotary inertia 
or curvature. This equation, obtained from reference 12 is 


(t) - [(iTb) 2 * -’ft * T^r) 


where x is the ratio of the sum of the face bending stiffnesses to the panel 
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T 


bending stiffness, 
becomes simply 


For identical faces of thickness tf = tf^ = tf^. 


T 


1 ( t f /h c \ 2 
3 \1 + t f /h c / 


(3a) 


The quantity £ is due to 
characterized by the ratio 


the orthotropic shear flexibilities of the core, 
Y = G c ^/G C y, and is given by 



(3b) 


The third frequency equation accounts for both the face bending stiffness 
and panel curvature but does not incorporate the effects of rotary inertia or 
orthotropic shear moduli of the core. This equation, obtained from the curved 
panel buckling solution of reference 10, is 


(£) 


2 - [(a7b) 2 



C 2 [m/ (a/b) ] 


1 + r{[m/(a/b)] 2 + n 2 }/ {[m/(a/b)] 2 + n 2 } 2 


( 4 ) 


Note: The last term in equations (11) and (13) of reference 10 contains a 

misprint, which has been corrected in equation (4) 2 above. 

It should be noted that equations (1), (2), and (4) are consistent. 
Equations (1) and (2) are identical when the effects of rotary inertia, curva- 
ture, and face bending stiffness are simultaneously neglected (X = C = x = 0) . 
Equations (1) and (4) are identical in the case of an isotropic core 
(r x = ry = r) when rotary inertia and the face bending stiffness are simulta- 
neously neglected (X = x = 0) . Equations (2) and (4) are identical in the 
case of an isotropic core (r v = r v = r, y = 1) when the panel curvature is 
neglected (C = 0). 7 

The basic assumptions made in the theories on which equations (1), (2), 
and (4) are based are: 

1. Core and faces show linear, elastic behavior. 

2. The transverse deflection of the panel comprises deformations due to 
both transverse shearing forces and bending moments. 


2 The symbols C, r, and x in this report are called Z^/tt 2 , , and S, 

respectively, in reference 10. 
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3. The bending stiffness of the core is negligible compared to the over- 
all bending stiffness of the panel, allowing the panel's flexural deformations 
to be described by an isotropic bending stiffness D s . 

4. Axial and circumferential inertia loadings are neglected. 

5. Transverse shear strains in the faces and all normal strains in the 
z direction are neglected. 

6. The cores are homogeneous; thus, for cellular type cores such as 
honeycomb the wavelengths of vibration must be at least several times the cell 
size . 


7. For the cylindrical ly curved panel, three additional assumptions are: 

a. The total panel thickness is small compared to the radius of 
curvature . 

b. The effect of the transverse shearing force Qy is neglected in 
the equation of circumferential equilibrium. 

c. The moment-distortion (curvature) relations are the same as for 
the flat panel. 

Although the natural frequencies predicted by equations (1) , (2) , and (4) 
are valid only for simply supported panels, the equations for modal density are 
applicable to panels having other boundary conditions: Once the frequency 

range of the lowest several modes is exceeded, the boundary conditions gener- 
ally have a negligible effect on the number of modes existing in a given fre- 
quency interval. This has been observed experimentally for flat panels and 
cylindrical shells that were essentially rigid in shear (refs. 13, 14, and 15). 
For mathematical arguments, see references 16 and 17. 

Modal Density Formulation 

For a given panel configuration (a/b, C, X, r x , ry, x, y, and to 0 fixed) 
equation (1) , (2) , or (4) represents a set of constant frequency curves in the 
m/ (a/b) , n plane. Figure 2 is a typical curve. Each combination of m and n 
represents the only mode that occurs in an area l[l/(a/b)] in the first quad- 
rant of the m/ (a/b) , n plane. The total area occupied by the N modes 
bounded by the curve to/to 0 = constant is approximately N[l/(a/b)]. This area 

, r ® 2 9 

is also given, approximately, by (1/2) J p z (w,0)d0, where 6 is defined in 

Jqi 

figure 2 and p 2 = [m/(a/b)] 2 + n 2 is determined in terms of to and 6 from 
equation (1), (2), or (4). The limits of integration are such that p 2 is 
real and positive in the region 0 <_ 0 <_ (tt/2) . As in reference 16, the approx- 
imate number of modes existing below the specified frequency to = constant is 
then expected to be asymptotically equal to 

N - \ £ f d2 P 2 (co,0)d0 

• / 01 
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The average number of modes AN existing in a frequency band Aw about 
the center frequency w defines the average modal density AN/ Aw. This, then, 
is approximated by 


AN. dN, 1 a 
Aw dw 2 b 



3p 2 (w,8) 
3w 


d0 + p 2 (w,02) 



dOj" 

p2(w ’ 6i:) d^T 


where Leibnitz's rule has been used. For the panel configurations considered 
herein the last two terms in the above equation turn out to be zero. 

If classical plate theory is used to predict the modal density of a flat 
sandwich panel, the frequency independent value dN/dw - (tt/4) (a/b) (l/w Q ) is 
obtained (ref. 17). On the basis of this theory, the modal density is seen to 
depend only on the panel 's lat eral dimensions, bending stiffness, and mass per 
unit area (w Q = (7r 2 /b 2 ) /D s /M) . For comparison with this classical plate 
theory value, the sandwich panel modal density estimates presented herein are 
all expressed in the dimensionless form 


1 dN 

(tt/4) (a/b) d(w/w Q ) 



9p 2 (w/w q ,0) 
3 (w/ w D ) 


(5) 


From the results that follow it can be noted that the right-hand side of 
equation (5) is independent of the panel length and width. These lateral 
dimensions appear only on the left-hand side of the equation, and since the 
term (a/b) (l/w Q ) is proportional to the product ab, the modal density AN/Aw 
is directly proportional to the surface area of the panel. Although equa- 
tions (1), (2), and (4) were derived on the basis of a rectangular planform, 
the modal densities obtained from them should be applicable to panels with 
other shapes of surface area S by replacing the rectangular area ab with 
the area S (ref. 17). 


ANALYSIS 


To calculate the modal density from equation (5) requires that 
3p 2 /3(w/w Q ) be expressed analytically as a function of w and 0. Unfortu- 
nately, p 2 appears as a cubic term in equations (1) and (4) and as a quartic 
term in equation (2), and thus is not readily obtained analytically. However, 
for the isotropic core (r x = ry = r) equations (1) and (2) can be solved 
exactly for 3p 2 /3(w/w Q ). Also, for an isotropic core, equation (4) can be 
used to obtain an approximate solution giving the combined effect of shear 
flexibility, panel curvature, and face bending stiffness. This is done in the 
following sections. It is then shown, by neglecting the face bending stiff- 
ness, that the isotropic results can be applied to panels with moderately 
orthotropic cores if "effective" shear flexibility and curvature parameters 
are used. Finally, an approximate solution is obtained for a flat panel, 
which gives the combined effect of orthotropic shear moduli and face bending 
stiffness . 
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Flat Sandwich Panel With Isotropic Core 


Face bending stiffness neglected- The frequency equation for this case is 
obtained from equation (1) by setting C = 0 and r x = ry = r, which results in 
two separate factors either of which is a solution if set equal to zero. 



1 + r 


[.■Hi 


•©']] - 4 



(6a) 


(6b) 


These equations are equivalent to those obtained in reference 18 and for a 
given value of p 2 yield three values of m 2 . The vibrations corresponding 
to the two higher frequencies, given by equation (6a) and the larger frequency 
root of equation (6b), are described in reference 18 as the thickness-twist 
and thickness-shear modes, respectively. The lower frequency vibration, given 
by the smaller frequency root of equation (6b), corresponds to a primarily 
flexural type of motion. This "bending set" of modes usually is of most 
interest in structural response problems and is described by 


20 2 - (r * x) fe) 2 



O 2 ] 2 * <t) 2 [‘ - 


(7) 


Since p 2 is independent of 0 (i.e., the curve w = constant is a circle in 
the m/(a/b), n plane), the limits of integration in equation (5) are 0j = 0 
and 02 = tt/ 2. Integration yields the following estimate for the modal density 
of a flat sandwich panel having an isotropic core: 


1 dN 

(tt/ 4) (a/b) d(a)/oj 0 ) 



1 + (1/2) [1 - (X/r)] 2 [r(a)/u 0 )] 2 
/l + (1/4) [1 - (X/r) ] 2 [r (w/ u> 0 ) J 2 


( 8 ) 

The variation in modal density with r(u)/w 0 ), as given by equation (8), 
is shown in figure 3 for X/r = 0 (i.e., rotary inertia neglected) and 
X/r = 0.279. For r = 0 (infinite shear stiffness) the frequency- independent , 
classical plate theory value for modal density is obtained. When shear flex- 
ibility is not neglected (r > 0) , the modal density increases with increasing 
frequency and with increasing shear flexibility. The value x/r = 0.279 cor- 
responds to the extreme case of a solid homogeneous panel whereas X/r <0.1 
is more typical of sandwich panel construction. Since the two curves in fig- 
ure 3 differ at most by about 6 percent, it is concluded that for the small 
values of X/r usually encountered in sandwich panels, the effect of rotary 
inertia on the modal density is negligible in comparison to the effect of shear 
flexibility. 
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The curves given by equation ( 8 ) are asymptotic to the curve 


1 dL = 2 (r-^l f 9 i 

0/4) (a/b) d(w/w 0 ) \ 0 ) o / J 

and for r(m/a) 0 ) >_ 3 .this straight line approximates the solid curve of fig- 
ure 3 to within 1 percent. However, equation (9) implies that the modal den- 
sity AN/A a) becomes independent of the face material properties at large 
values of r(u)/w 0 ). This physically unrealistic result arises because equa- 
tion ( 8 ) does not account for the bending stiffness of the panel faces. 

Effect of face bending stif fness- In view of the results of the previous 
section it is assumed here that rotary inertia can be neglected. Equation (2) 
can then be used to determine the effect of face bending stiffness on the 
modal density. For an isotropic core, equation (2) can be written as 


r 


w 



rp 



1 

1 + rp 2 


which is independent of 0 . Use of equation (5) yields 


( 10 ) 


1 dN a 1 

Ii d /_£_\ 3 (m/mq) 

4b \“o/ 3p 2 


/l + rp 2 


J\ + t( 1 + rp 2 ) 


1 - 


rp 


2(1 + rp 2 ) [1 + x (1 + rp 2 )] 


( 11 ) 


for the estimate of modal density. The variation in modal density with 
r(o)/co 0 ) is obtained from equations (10) and (11) by treating t as a param- 
eter and the quantity rp 2 as a variable. For the case of identical face 
sheets, this variation is shown in figure 4 for various values of the face to 
core thickness ratio tf/h c . It is seen that neglecting the face bending 
stiffness (tf/h c = 0 curve) leads to an overestimate of the modal density by 
an amount that increases with increasing frequency and with increasing ratios 
of face to core thickness. However, for many sandwich configurations tf/h c 
is on the order of 0.1 or less and often the quantity r/co 0 is so small that 
a frequency of several thousand Hertz is required to produce a value of 
r(w/w 0 ) on the order of 1.0. Thus, there are many practical situations where 
the face bending stiffness can be neglected when computing the modal density. 

When the face bending stiffness should be taken into account, the follow- 
ing simple equation is useful 


8 


dN 


( 12 ) 


2r(o)/a> 0 ) 


(tt/4) (a/b) d(u,/u> 0 ) “/i 7 T [2r((0/aSo)] : 


For r(w/u) 0 ) _> 3 and tf/h c £ 0.'2 (identical face sheets), this result predicts 
modal density values that are within about 1 percent of the values given by 
the exact solution. (Eq. (12) is obtained from eqs . (10) and (11) by taking 
rp 2 >> 1.) 


Flat Sandwich Beam 


The differential equations in appendix B can be applied to a flat beam 
whose length b is in the y direction if all terms that involve the x 
coordinate are dropped. This has the effect of eliminating the thickness- 
twist solution and all terms involving m/(a/b) from the bending and thickness- 
shear solutions. Thus, the frequency equation for the beam "bending" modes 
can be obtained directly from the panel "bending" solution by setting the 
m/ (a/b) terms equal to zero. The modal density is then simply obtained from 
dN/d(m/o) 0 ) = dn/d(to/m 0 ). 


Face bending stiffness neglected - This case is obtained by setting 
m/(a/b) = 0 in equation (7). Differentiating with respect to w/m 0 yields 




dN 


* » ^ t) 2 

(> - f )> 




(13) 


The variation in modal density with frequency, as given by equation (13), 
is shown in figure 5 for various values of the shear flexibility parameter r 
and for two values of x/r. As in the case of panels, realistic values of 
X/r give essentially the same results as X = 0 (rotary inertia neglected). 
The r = X = 0 curve corresponds to classical beam theory, which predicts that 
modal density decreases with increasing frequency. When shear flexibility is 
accounted for (r > 0) , equation (13) predicts that the beam modal density is 
nearly constant over a wide frequency range. However, this result is valid 
only for conditions where the face bending stiffness is unimportant. 

Effect of face bending stiffness- If rotary inertia is neglected, 
equation (2) can be used to determine the effect of the face bending stiffness 
on a beam's modal density by letting m/(a/b) = 0. The result is 


0) 




1 

1 + rn 2 


(14) 
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dN 


/l + m 2 


(15) 


2(1 + rn 2 ) [1 + t (1 + rn 2 )] ) 

The variation in modal density with frequency as given by equations (14) 
and (15) is shown in figure 6 for several values of r and with a face-to- 
core thickness ratio of 0.05 (identical face sheets). Comparison of the 
curves in figure 6 with the corresponding curves in figure 5 reveals that neg- 
lect of the face bending stiffness again leads to an overestimate of the 
modal density by an amount that increases with increasing frequency. The 
face bending stiffness is also seen to have a smaller effect on beams with 
stiff cores (small r) than on beams with more flexible cores (larger r) . 

Figures 5 and 6 both indicate a large rise in modal density as w/u) 0 
approaches zero. However, it must be remembered that the equations presented 
for modal density are estimates that have been obtained by representing the 
number of modes by a continuous function of frequency. Such a representation 
is not really meaningful at frequencies near that of the fundamental mode 
o)/o) 0 = 1 as there are usually too few modes involved. This is illustrated in 
figures 5 and 6 by the line labeled N 6. For frequencies to the left of 
this line only five modes exist. Thus, the large rise in modal density as 
m/u o -* 0 is physically unrealistic. A similar situation occurs in the modal 
density estimates for curved panels. 


(£) 


d(u/u 0 ) 


dn 


2n 


/T 


t( 1 + rn 2 ) 


Curved Sandwich Panel With Isotropic Core 

Face bending stiffn ess neglected- When the parameter C is retained in 
equation (1), the effects of panel curvature can be determined. For an iso- 
tropic core (r = r x = ry) , the "bending set" of modes is described by 

2p 2 = (r + X)^— ^ - rC 2 cos 4 (6) 

V^o/ 


+ 



rC 2 





The modal density is then given by 


(16) 


1 

— fL 
4 b 






»(■- DO -Dfr ff 



( 17 ) 
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II 


where 


a = 1 



cos 4 (0) 


The lower limit of integration is 


if 


and 






1 


are both satisfied; otherwise 0^ = 0. 

When X = r = 0 (rotary inertia and shear flexibility neglected) , equa- 
tion (17) can be expressed in terms of the complete elliptic integral of the 
first kind (ref. 16). For nonzero values of X and r, equation (17) is 
expressible as a hyperelliptic integral. Some hyperelliptic integrals can be 
reduced to the sum of elliptic integrals; however, the results presented 
herein were obtained by numerically integrating equation (17). These results 
are presented in figure 7, where the variation in modal density with (uj/u 0 )/C 
is shown for various values of the product rC. The numerical results shown in 
figure 7 are all for X = 0 (rotary inertia neglected). 3 However, in view of 
the relatively insignificant effect of rotary inertia on the modal densities 
of the beam and flat panel, it is likely that the curves in figure 7 would be 
only slightly altered by realistic values of X/r. 

The peaks in modal density about uj/u) 0 = C (the so-called "ring fre- 
quency") are due to a concentration of modes occurring near this frequency and 
is most pronounced for the r = 0 case. However, the sharp spike is due to 
the singularity in equation (17) at u)/u) 0 = C (when X = 0) and is physically 


3 For X = 0 and (u)/w 0 )/C <_ 1, the integrand of equation (17) diverges at 
01, but for (w/w Q )/C ^ 1, the integral is convergent. In the numerical inte- 
gration of equation (17) the contribution of the integrand in the neighborhood 
of 0} was approximated by replacing the integrand with the first three terms 
of its Taylor series representation. For r = X = 0, the results obtained in 
this manner agree within 0.5 percent of the results given by the elliptic 
integral solution. 
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unrealistic. A better estimate for the peaks in modal density about the ring 
frequency can be obtained by directly calculating the number of curved panel 
modes ANq occurring there from the frequency equation 



(Eq. (18) is obtained from eq. (16) or (4).) 

For r = 0, the results of such calculations for the interval 
0.95C < w/u) 0 < 1.05C are presented in table 1 for various degrees of curvature 
(a/b was taken equal to 4 / tt) . The number of modes ANp that the correspond- 
ing flat panel has in the same frequency interval is also shown for comparison 


TABLE 1.- COMPARISON OF CURVED AND FLAT PANEL MODAL DENSITIES AT THE 

RING FREQUENCY (r = X = 0) 






Number of modes in 


Curvature 

parameter 

Frequency interval 
n_Q5n < < i.05r 

frequenc 

y interval 

[> 

n 

C 



o 

Flat panel 

Curved panel 

ANp 





ANp 

AN C 


50 

47.5 


52.5 

5 

5 

1.0 

75 

71.25 

-y 

78.75 

6 

8 

1.33 

100 

95 

-y 

105 

10 

13 

1.30 

200 

190 

-y 

210 

21 

27 

1.29 

300 

285 

-y 

315 

29 

43 

1.48 

400 

380 

-y 

420 

39 

55 

1.41 

500 

475 

-y 

525 

49 

76 

1.55 

750 

712.5 

-y 

787.5 

76 

115 

1.51 

1000 

950 

- y 

1050 

101 

155 

1.53 


For 300 <L C _S 1000, the actual peak in curved panel modal density is about 
1-1/2 times that of the flat panel value. This result is indicated in fig- 
ure 7. For 75 <_ C < 300, a somewhat smaller increase is evident. For C 
less than about 50 the effect of curvature is so small that the panel essen- 
tially behaves as if it were flat, not only at 0.95 < (a)/u 0 )/C < 1.05 but over 
the entire frequency range. 

Except in the r = 0 case, figure 7 obscures the fact that when u)/w 0 
exceeds C the curved panel modal density becomes asymptotic to that of the 
corresponding flat panel (C = 0) . This behavior is more clearly seen by 
replotting figure 7 in terms of r(u)/u) 0 ) for fixed values of rC as shown by 
the solid curves in figure 8. The curve labeled rC = 0 denotes the flat 
panel solution. Because of the singularity in equation (17), the modal den- 
sities at u)/oj = C were estimated from equation (18) in the same manner as 
was done for r = 0. There is essentially no difference, at or above the ring 
frequency, between the curved and flat panel results whenever the product rC 
is greater than two. 
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Useful approximations to the solid curves in figure 8 are given by the 
equation 


dN 


sin -1 ^/r(m/ aj Q )/rC 


Or/4) (a/b) d(aj/a» 0 ) 


= ^ 


2r — 


2r JiL 


tt/2 


0) r 


for 


for 


w/o) 




< 1 


> 1 


> (19) 


Results from equation (19) are shown by the dashed curves in figure 8. It is 
seen that this simple approximation to the solution of equation (17) is fairly 
accurate for values of rC greater than about 3 if r(u>/u) 0 ) is greater than 
about rC/3. (Eq. (19) is obtained by taking rp 2 >> 1 and X = 0 in eq. (16) 

with Gj = cos -1 / (m/( 0 o )/C and 02 = tt/ 2 in eq. (5).) 


Effect of face bending stiffness- Since equation (17) does not account 
for the face bending stiffness, it undoubtedly overestimates the modal density 
at the larger values of r(aj/u 0 ). (The curves in fig. 7 correspond to 
0 <_ r(w/u) 0 ) <.10.) At the larger values of r(u)/u> 0 ) an approximate solution 
giving the combined effect of panel curvature and face bending stiffness can 
be obtained from equation (4) by neglecting 1 with respect to rp 2 as -in the 
derivation of equations (12) and (19) . This leads to 


dN 


2r(tu/w 0 ) 


(tt/ 4) (a/b) d( u/Ub ) J Y + T[2r(aj/a)o)] 2y e ( , /2) /I 


I 


TT/2 


d0 


6 2 cos 4 0 


(20a) 


where 


and 


t ( 2rC) 2 

1 + x [2r (w/ w 0 ) ] 2 



if 



if 



(20b) 


y (20c) 


By making the substitution t 
of equation 254.00 of reference 19 


= cos 2 0, equation (20a) takes the form 
and can be expressed as follows. 
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For 6 < 1: 


where 


and 


where 


and 


dN 


2r{>/io 0 ) 


F (<h >kj) 


0/4) (a/b) d(«,/u> 0 ) „ t [ 2r (oj/too) ] 2 "' 2 /TVJ 

. / 26 
kl_ /l + 6 


if 


C - 


> 1 




sm 


_j /[(aj/a> 0 )/C](l + 6) 


1 + 6 [ (w/w 0 )/C] if 


U)/u) 


< 1 


For 6 1 : 

1 


dN 


2r(w/w 0 ) F(<J> 2 ,k 2 ) 


(t/ 4) (a/b) d(«/%) /H. t Ur (»/.*>)] 2 


it/ 2 


/2S 


a + 6 


26 


<j >2 = sin 


-1 


26 [ (u)/a) 0 )/C] 


1 + 6[0/co 0 )/C] 


> (21a) 


j 


L (21b) 


J 


In equations (21a) and (21b), F(<J>,k) is the incomplete elliptic integral of 
the first kind. 


The variation in modal density with r(w/to 0 ) as given by equations (21) 
is shown by the solid curves in figures 9(a) through 9(c) for values of rC 
equal to 3, 5, and 10, respectively. Each figure is for identical face sheets 
with a face-to-core thickness ratio of tf/h c = 0.1. 

The exact solution for the modal density based on the entirety of equa- 
tion (4) has not been obtained. Without having this exact solution, the 
errors introduced by the approximation rp 2 >> 1, on which equations (21) are 
based, cannot be firmly established. However, for C = 0 (flat panel) and 
x = 0 (face bending stiffness neglected), equations (21) reduce to the approx- 
imate solutions given by equations (12) and (19), respectively. The latter 
two approximate solutions were shown to agree well with the corresponding 
exact solutions (eqs. (11) and (17), respectively) over certain parameter 
ranges. This suggests that the solution given by equations (21) may also be 
a fairly reliable estimate of the modal density for the same range of 
parameters . 
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In the case of the flat panel, the effect of the face bending stiffness 
is relatively small for r(co/a) 0 ) < 5 and tf/h c <. 0.1 (see fig. 4). Thus, if 
the ring frequency is in this range (rC < 5), the modal density below the ring 
frequency would be expected to be significantly affected by the curvature and 
only slightly affected by the face bending stiffness. The approximate solu- 
tion given by equations (21) does predict this behavior as can be seen from 
figures 9(a) and 9(b) (rC = 3 and 5, respectively). Below the ring fre- 
quency the approximate solution (solid curve) closely follows the exact solu- 
tion (dot-dashed curve) for the curved panel without face bending stiffness 
given by equation (17). Above the ring frequency, the approximate solution 
closely follows the exact flat panel solution (dashed curve) given by 
equation (11). 

Next, consider the ring frequency in the frequency range where the face 
bending stiffness of the flat panel becomes important (rC > 5 for t£/h c = 0.1). 
In this situation, the modal density below the ring frequency would be 
expected to increase at a less rapid rate than given by equation (17) . This 
too is predicted by equations (21) as shown in figure 9(c) (rC = 10) by the 
increased separation between the solid and dot-dashed curves in the range 
5 <. r(w/u) 0 ) < 10. Above the ring frequency (r(w/w 0 ) > 10 in fig. 9(c)), the 
approximate solution again approaches the exact flat panel solution. 

As a further check on the solution given by equations (21) the modal 
density for rC = 10 was determined directly from equation (4) by counting the 
actual number of modes occurring in intervals of A[r(w/to 0 )] = 1.0. The 
values for the modal density obtained in this manner are shown by the circular 
symbols in figure 9(c) and are quite close to the estimate (solid curve) given 
by equations (21) . (The parameter values used for this computation were 
a/b = 1.27, C = 200, r = 0.05, and t f /h c = 0.1.) 


Effect of Core Orthotropy 


Face bending stiff ness neglected- When the effects of face bending stiff- 
ness and rotary inertia are neglected, the equation governing p 2 for a 
curved sandwich panel having an orthotropic core is obtained from equation (1) 
by setting X = 0. This yields 



C 2 cos ^ ( 0 ) 



C 2 


cos 4 (0) 


(‘ 



cos 2 (0) 


0 ( 22 ) 


Except when y = 1 (isotropic core), equation (22) is cubic in p 2 and 
no attempt is made to solve it in this form. However, an indication of the 
effect of y can be obtained by considering the conditions 


lyP' 


» 1 


r yP‘ 


» 1 


(23) 
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Under these conditions P 2 is approximately equal to 


r [0/to 0 ) 2 - C 2 cos 4 (0)] 

— ^ > 0 (24) 

y(l - [ Cy - l)/y]sin 2 (0) } 


Substituting into equation (5), and noting that ©i = 0 and 02 = ff/2, yields 


1 dN a 7 / V\ / J?_\ 

Ctt/ 4) Ca/bD d (to/ w 0 ) “ 


(25) 


Comparison of this result with equation (9) suggests that the isotropic 
results can perhaps be applied to panels with orthotropic cores if r is 
replaced with an effective shear-flexibility parameter r e ££ = r y / /f . However, 
it is also likely that C will have to be replaced with an effective curva- 
ture parameter since equation (25) is not valid where the effects of curvature 
are expected to be most pronounced (to/u> 0 <_C). 4 

This idea was examined by solving equation (22) for w/w 0 in terms of m 
and n: 


/ m\ 2 _ { [m/ (a/b) ] 2 + n 2 } 2 + C 2 [m/(a/b)] 4 

1 + ? {[m/(a/b)] 2 + n 2 } 2 


where ? is the same as in equation (3b). 

For specified values of a/b, r y , y, p, and C, the number of modes N 
existing below any frequency «/u 0 can be calculated from equation (26) and 
a plot of N vs. w/u) 0 is readily constructed as shown in figure 10. The pos- 
sibility of duplicating this plot (hence the modal density) by calculating N 
from the isotropic frequency equation (eq. (18)) with r replaced by some 
r eff an< ^ with C replaced by some C e ff was then considered; that is, from 


_w\ 2 = { [m/ (a/b) ] 2 + n 2 } 2 + C eff [m/ (a/b) ] 4 

1 + r eff { [m/(a/b)] 2 + n 2 } {[m/(a/b)] 2 + n 2 } 2 


(27) 


The two curves in figure 10 were obtained from the exact equation 
(eq. (26)) and show that, for r y //y~ held constant, a decrease in y (stif- 
fening of core in circumferential direction) causes a decrease in the cumula- 
tive number of modes. This same effect is produced by eqaution (27), for 
r e ff held constant, if C e ff is made to increase as y decreases. For 


4 The inequality p 2 > 0 expressed by equation (24) cannot be satisfied at 
0 = 01 if w/u) 0 <_ C. 
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values of y from 0.4 to 2.5 (representative of honeycomb cores) with 
1 <_ a/b <_ 2.55, r y //7 up to 0.125, and C up to 500, it was found that the 
plots of N vs. aj/m 0 obtained from equation (26) were duplicated quite well 
by equation (27) if r g ££ and C e ££ were taken as 


_ r y _ *1 

r eff " Sy " b 2 


y°Qx D Qy 


(28a) 


c eff “ < 




( y ) 0 - 15 

(C 


if 


if 


C - 
U)/ id 


< 0.9 


(28b) 


> 0.9 


If the above expressions for r e ££ and C e ££ are used, equation (27) 
(empirical) predicts modal densities that are within 5 percent of the results 
predicted by equation (26) (exact) for a flat panel and also for a curved 
panel when (id/<d 0 )/C is greater than 0.9. In the range 0.1 <_ (u)/o> 0 )/C <_ 0.90, 
the accuracy of equation (27) depends on r g ££C. For r e ffC <_ 2, the average 
error in the modal density is less than 5 percent; as r e ffC increases to 5, 
the average error increases to about 10 percent, and for r e ffC = 10, the 
average error is about 15 percent. Typical results are shown in figures 11(a) 
and 11(b). Figure 11(a) shows exact plots of N vs. co/ui 0 for C = 0 and 
C = 100 with Ty = 0.00707, y = 1/2, and a/b = 2.55. The corresponding plots 
obtained from equation (27) with r e ££ and C e ff defined by equations (28) are 
shown in figure 11(b). 


The form of equation (27) is exactly the same as equation (18). Thus, 
when rotary inertia and face bending stiffness effects are small, the modal 
density estimates obtained for panels having isotropic cores appear applicable 
to panels having moderately orthotropic cores such as honeycomb by simply 
replacing r and C with the quantities r g ££ an< ^ C e ££, respectively, defined 
by equations (28) . Note that the "effective" shear stiffness is the geometric 
mean of the two orthotropic stiffnesses. 

Another check on the validity of this empirical approach is that for a 
flat panel, equations (27) and (28) predict the variation in N with w/w 0 
to be unaffected by a 90° core rotation. This has been verified by using 
equation (26) (the exact equation) to show that, for the range 0.4 <_ y <_ 2.5, 
plots of N vs. id/idq are indeed unaffected by such a rotation. A typical 
example is shown in figures 12(a) and 12(b). Figure 12(a) is for a flat panel 
(C = 0) with Ty = 0.00707, y = 1/2, and a/b = 1.27. Figure 12(b) corresponds 
to the same flat panel but with the core rotated 90° (r y = 0.01414, y = 2) . 

The resulting curves of N vs. id/id q are virtually identical. 

Effect of face bending stiffness- In the case of a flat panel an approxi- 
mate solution giving the combined effect of orthotropic core shear moduli and 
face bending stiffness on the modal density can be obtained by neglecting 
certain terms in equation (2) . If the conditions given by equations (23) are 
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satisfied, then by order-of-magnitude considerations, equation (2) simplifies 
to the approximation 

r y Tp 4 + [1 + (y - l)cos 2 0]p 2 - r y(^) = 0 (29) 

Solving for 9p 2 /9 (a)/u) 0 ) , substituting into equation (5), and noting that 
0! = 0 and 0 2 = tt/2, yields 


1 dN 

(tt,/ 4) (a/b) d(a)/w 0 ) 


2r y (m/o) o ) 

tt/2 



TT/2 


d0 


/ T [2r y (w/ oj 0 ) ] 2 + [1 + 


(y - l)cos 2 


0]2 

(30) 


Making the substitution t = - (y - l)cos 2 0 reduces equation (30) to 
the form of equation 259.00 of reference 19 so that the solution is 


where 


(ir/4) (a/b) d ( oj/oj, 


* ( 2r ««^) el 


/4 K(k) 
tt/ 2 


eff 


y _ ir 

^ ' b 3 


/ d Qx d Q> 


2k 2 


[l ♦ ^ (2r eff Jl-) 2 ] [ T * T ( 2 r eff ^) 2 ] 
‘ 1 ' [ ltT ( 2r ^)] /f 


and K(k) is the complete elliptic integral of the first kind. 






(31) 


V (32) 


The terms in equation (31) that involve the orthotropic shear moduli are 
3 and ? e f£- Both of these terms, and therefore the modal density given by 
equation (31) , are unaffected by a 90° rotation of the core . This is shown in 
figure 13 by the dashed curves labeled y = 1/2 or 2, y = 1/4 or 4, and 
y = 1/8 or 8. The solid curve is for an isotropic core (y = 1) and is 
obtained from the exact solution (eq. (11)). All the curves shown in fig- 
ure 13 are for identical faces with a face-to-core thickness ratio of 


tf/h c = 0.1. 


Note that equation (31) depends on both J Dq^Dq^ and y = Dq /Dq . 

It is only when the face bending stiffness is negligible that the effect of 
an orthotropic core on the modal density of a flat panel ca n be desc ribed 


solely in terms of the effective shear stiffness D, 


Qeff 


/ 'VQj 
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For y = 1 (isotropic core), equation (31) reduces to equation (12). 

This suggests that for the moderately small range of y usually encountered 
in honeycomb cores (0.4 < y < 2.5), equation (31) should be fairly accurate if 
tf/h c <.0.2 and r e ff(w/co 0 ) >_ 3 are both satisfied. 


EXPERIMENT 


Few experimental data are available in the literature on the mode shapes 
and natural frequencies of vibrating sandwich panels, especially for frequen- 
cies substantially above that of the fundamental mode. Although about 35 
modes were excited in each of the two panels tested in reference 12, the panel 
configurations and frequency range were such that the maximum value obtained 
for r e ff(to/io 0 ) was about 0.25. At this small value of r e f p (w/w 0 ) , the shear 
flexibility of the core theoretically has a relatively small effect on the 
modal density (see fig. 3). Thus, it was necessary to conduct some experi- 
ments that would cover a higher range in r e ff(m/w 0 ) for the purpose of 
obtaining results that could be compared with theory. 


Apparatus 


Sandwich panels- The specimens tested consisted of four flat, rectangular 
sandwich panels constructed from aluminum honeycomb cores bonded to stainless 
steel face sheets with a structural adhesive. The effective length and width 
of each panel (measured between the boundaries of the panel support fixture) 
were a = 28.5 in. (72.4 cm) and b = 24.0 in. (61.0 cm), respectively. 

Table 2 lists the pertinent core and face properties along with the resulting 
panel parameters. The weight of the bonding material was lumped with the 
weight of the faces to produce an effective face density called (Pf) e ff 


Values of the core shear moduli were calculated from the following 
equations (ref. 20) 

G ""N 

_ _s_ [1 + (u/v) cos n] 2 

Lp Ps [1 + (U/v)] 2 

r 1 + (u/v)cos n L 

l D p s Pc 1 + (u/v) | 


r _ G s , (u/v) sin 2 n 

T p s p c i + (u/v) 


(33) 


The quantities G s = 3.85xl0 6 lb/in. 2 (26.5 GN/m 2 ) and p s = 2.50xl0 -i+ 
lb-sec 2 /in. 4 (2.67 Mg/m 3 ) are the shear modulus and density of the honeycomb 
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TABLE 2.- DESCRIPTION OF TEST PANELS 


Length: a = 28.5 in. (72.4 cm), each panel 
Width: b = 24.0 in. (61.0 cm), each panel 

Materials 

Panels 1 and 2: Faces of 302 stainless steel sheet 

Cores of 1/8-5052 - 0.002 aluminum honeycomb 

Panels 3 and 4: Faces of 321 stainless steel sheet 

Cores of 0X-3/ 16-5052 - 0.0007P overexpanded aluminum honeycomb 



Core properties 

Face properties 

Panel 


G 

cx 

G 

. . cy 

h 

c 

(p f } eff 

E 

f 

V 


m 


B 




in. 


Bflj 







1 

1.21xl0' 5 

129 

131xl0 3 

903 

S3.0xl0 3 

365 

0.502 

1.28 

80. 3xl0 -5 

8.58 

26x10® 

0.18 

0.245 

0.0154 

0.391 

2 

1.21xl0~ 5 

129 

53.0xI0 3 

365 

131xl0 3 

903 

0.502 

1.28 

81.4xl0' 5 

8.70 

26x10® 

0.18 

0.245 

0.0154 

0.391 

3 

0. 316 x 10 -5 

33.8 

21.7xl0 3 

150 

23.9xl0 3 

165 

0.503 

1.28 

80.3xl0' 5 

8.58 

28x10® 

0.19 

0.245 

0.0203 

0.516 

4 

0.316xl0‘ 5 

33.8 

23.9xl0 3 

165 

21 . 7xl0 3 

150 

0.503 

1.28 

81.1xl0 -5 

8.67 

28xl0 6 

0.19 

0.245 
1 

0.0203 

0.516 




Panel 

a/b 

t f/ h c 




Y 

1 

1.19 

0.0307 

738 

0.0221 

0.0449 

2.47 

2 

1.19 

0.0307 

734 

0.0221 

0.0450 

0.405 

3 

1.19 

0.0404 

843 

0.114 

0.00993 

0.9081 

4 

1.19 

0.0404 

840 

0.114 

0.00993 

1.10 

















4 



T 


L 


Sketch (a) Geometry of 
honeycomb cells. 


material , respectively; u, v, and n are 
defined in sketch (a) . The honeycomb shear 
modulus in the longitudinal (L) direction of 
the core is designated Gl and the shear 
modulus in the transverse (T) direction of 
the core is designated Gp. The first and 
second of equations (33) are lower and upper 
bounds, respectively, to the longitudinal 
modulus and the arithmetic average of these 
values was used for G^. 


Experimental values of shear moduli obtained from single-block shear 
tests of the overexpanded cores used in panels 3 and 4 were supplied by the 
honeycomb manufacturer. However, these tests showed considerable scatter with 
values of /G^G-p from 21 to 31 percent lower than the result obtained from 
equations (33). Since values of shear moduli obtained from block shear tests 
are known to be generally low in comparison with values obtained from theory 
and other test methods (ref. 21) , it is felt that the theoretical values used 
herein are more accurate than the experimentally obtained values. For the 
normally expanded cores used in panels 1 and 2, the handbook value of Jg lGt 
was within 2 percent of the theoretical result used. 


For the cores used in panels 1 and 2, n = 43° and u/v = 1.6. For the 
cores used in panels 3 and 4, n = 72° and u/v = 1.2. The approximate geometry 
of the honeycomb cores and the orientation with respect to each panel is shown 
in figure 14. Panels 1 and 2 were nearly identical except for a 90° rotation 
of the core. The primary difference between panels 3 and 4 was also a 90° 
core rotation. 



Sketch (b) Panel mounting fixture. 


Support and ex citation system- 
The panel mounting fixture consisted 
of two aluminum frames bolted to the 
outer 3 in. (7.62 cm) of the panel 
perimeter to provide a partially 
clamped edge condition (see 
sketch (b)). Each steel bolt was 
torqued to 175 in. -lb (19.8 m-N) . To 
prevent crushing of the core, the 
outer 3 in. (7.62 cm) of honeycomb 
cells were filled with a liquid 
aluminum potting material. 


-Panel 


1/32- in. (0.794 mm) 
thick rubber washer 


2 


To vacuum 
3/4- in. P um P 

(1.90 cm) dia. 


/777777777777ZW 


I Shaker motion 


Sketch (c) Vacuum cup attachment. 


The fixture was supported by a 
wooden frame, and the panel was 
excited from below by a small (25-lb 
force (111 N)) permanent magnet 
shaker. The shaker was coupled to 
the panel surface by means of a small 
vacuum cup to avoid the fastening of 
attachment points to the panel (see 
sketch (c)). This arrangement proved 
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to be quite satisfactory for transmitting motion to the panel and it allowed 
the excitation point to be easily changed. 

Instrumentation- The shaker was driven by the amplified signal of a 
variable frequency oscillator. A hand-held vibration pickup was used to 
detect the panel motion. The frequencies of the oscillator and vibration 
pickup signals were measured by an electronic frequency counter. These 
signals were also monitored on an oscilloscope in the form of Lissajous 
figures . 


Test Procedure 

For each panel, between 70 and 80 consecutive modes of vibration were 
excited by varying the frequency and location of the applied excitation. The 
mode shape at each resonance was visualized by the formation of Chladni fig- 
ures produced by the collection of sand particles (16 mesh size) along node 
lines. A few of these sand patterns are shown in figure 15. For some of the 
higher frequency resonances, the Chladni figures were not sharply defined. In 
these cases, the phase differences between the oscillator signal and the vibra- 
tion pickup signals were used to detect the node lines. In some instances, 
two or more modes occurred at nearly the same frequency, and careful position- 
ing of the excitation point was required to produce distinct Chladni figures. 


Test Results 

Frequencies are given in table 3 for the four test panels corresponding 
to the maximum resonance response for the modes listed. The mode number m 
in the x direction (n in the y direction) indicates m - 1 (n - 1) lines 
of zero deflection between the panel boundaries x = 0 (y = 0) and x = a 
(y = b) . For panels 1 and 2 the ratio Gp/G-j. was equal to about 2.5, which 
caused corresponding higher modes of these two panels to occur at signifi- 
cantly difference frequencies. For panels 3 and 4, Gp/G-p was nearly unity, 
and corresponding modes for these panels occurred at nearly the same 
frequency. 


COMPARISON OF THEORY AND EXPERIMENT 


The experimental values of { 1/ [ (ir/4) (a/b) ] } [AN/A (w/co 0 ) ] and their varia- 
tion with r e £f(w/w 0 ) are shown in figures 16(a) and 16(b) for panels 1 and 2, 
respectively, and in figures 16(c) and 16(d) for panels 3 and 4, respectively. 
The horizontal lines indicate the intervals of r e ff (w/a) 0 ) used to compute the 
corresponding data points. The solid curve in each figure is the theoretical 
estimate for modal density given by equation (8) with r replaced by r e £f. 
For all four panels, the quantity X/r e £f is much less than one; therefore, 
theoretically, the effect of rotary inertia is negligible. For the range of 
r eff(w/(jJ 0 ) and t£/h c covered by the experiments, the theoretical effect of 
the face bending stiffness is also negligible. 
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[Panels 


Panel 1 


Panel 2 


Panel 3 


Panel 4 


TABLE 3.- EXPERIMENTAL RESONANCE FREQUENCIES, H z 
1 and 2 heavy line encloses modes for which r e ff(w/(i) 0 ) < 1.1; panels 

3 and 4, r e ff(to/w 0 ) < 3.0] 


m/n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1 

294 

713 

1130 

1890 

2470 

3070 

3720 

4360 

5020 

5660 

6330 

2 

615 

950 

1425 

1960 

2590 

3230 

3850 

4460 

5140 

5770 


3 

1018 

1360 

1750 

2280 

2780 

3390 

3950 

4610 

5280 

5910 


4 

1600 

1840 

2200 

2620 

3180 

3790 

4310 

4920 

5520 



5 

2162 

2460 

2815 

3130 

3630 

4130 

4750 

5300 

5900 



6 

2980 

3170 

3440 

3750 

4200 

4720 

5220 

5740 




7 

3716 

3870 

4010 

4370 

4780 

5230 

5750 

6220 




8 

4470 

4660 

4930 

5090 

5460 

5890 






9 

5260 

5430 

5630 

5820 

6160 







10 

6070 











1 

300 

726 

1180 

2260 

3060 

3960 

4900 

5940 




2 

530 

991 

1560 

2350 

3170 

3990 

4970 





3 

940 

1300 

1780 

2450 

3270 

4120 

5080 





4 

1340 

1690 

2120 

2700 

3530 

4360 

5150 





5 

1870 

2020 

2520 

3080 

3790 

4570 

5370 





6 

2290 

2510 

2890 

3440 

4090 

4820 

5630 





7 

2750 

2930 

3360 

3840 

4470 

5150 

5980 





8 

3220 

3420 

3800 

4250 

4800 

5540 






9 

3710 

3960 

4240 

4610 

5220 

5870 






10 

4200 

4421 

4680 

5170 

5630 







11 

4610 

4910 

5070 

5640 

6070 







12 

5130 

5360 

5590 

6080 








13 

5700 

5800 

6100 









1 

280 

674 

940 

1360 

1760 

2140 

2520 

2910 

3290 

3640 


2 

443 

760 

1070 

1430 

1810 

2200 

2630 

3020 

3360 



3 

810 

990 

1300 

1570 

1930 

2310 

2700 

3080 

3430 



4 

1010 

1230 

1500 

1770 

2110 

2440 

2800 

3140 

3500 



5 

1330 

1530 

1750 

1980 

2310 

2580 

2930 

3280 




6 

1760 

1830 

2020 

2220 

2500 

2780 

3130 

3440 




7 

2000 

2140 

2290 

2490 

2700 

3000 

3350 





8 

2360 

2440 

2560 

2750 

2960 

3210 

3540 





9 

2650 

2710 

2890 

3060 

3230 

3540 






10 

2950 

3050 

3200 

3300 

3500 







11 

3300 

3390 

3460 









12 

3590 











1 

288 

638 

951 

1380 

1780 

2140 

2540 

2910 

3310 

3670 


2 

437 

754 

1060 

1430 

1820 

2200 

2620 

3060 

3380 



3 

784 

990 

1260 

1580 

1940 

2300 

2700 

3100 

3460 



4 

981 

1220 

1480 

1780 

2100 

2450 

2800 

3130 

3590 



5 

1310 

1540 

1720 

2010 

2310 

2630 

2940 

3310 




6 

1780 

1820 

2000 

2230 

2530 

2810 

3140 

3480 




7 

2030 

2130 

2290 

2490 

2720 

3010 

3360 





8 

2340 

2430 

2580 

2760 

2970 

3240 

3530 



- 


9 

2630 

2700 

2870 

3030 

3260 

3500 






10 

2930 

3050 

3210 

3340 

3540 

3730 






11 

3270 

3390 

3460 

3620 








12 

3650 

3700 
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Comparison of the results presented in figure 16 indicates that the 
theory gives a fairly good estimate of the average modal density except at the 
smaller values of r e ff (to/to 0 ) . The values of r e ff(a)/u) 0 ) at which the compar- 
ison becomes poor correspond to the frequency range in the vicinity of the 
fundamental mode where relatively few modes occur. In this frequency range, 
the continuous frequency representation of a discrete number of modes is 
unrealistic. For example, the first two test points in figure 16(c) represent 
a total of only five modes. In contrast, the last test point in figure 16(c) 

represents a total of 24 modes. (Note that this is six times the number of 

modes predicted by classical plate theory.) 

Panels 3 and 4 are nominally identical, the major difference being about 
a 20 percent variation in y. However, comparison of figures 16(c) and 16(d) 
shows that the experimental results obtained from panel 3 fall very close to 
the theoretical curve at nearly every point while the results from panel 4 are 

more scattered. This is due to a "clumping" of more than the average number 

of modes (as predicted by theory) in one frequency interval while an adjacent 
frequency interval has fewer than the average number of modes. (This clumping 
effect is also seen in figs. 12(a) and 12(b), for example, at 170 <u)/w 0 <180.) 
Depending on the size of frequency interval chosen, this clumping may or may 
not produce noticeable variations in modal density from the average. The 
theory gives only the average modal density and does not predict variations 
from this average. 

It should also be noted that the test panels were fastened in a semi- 
clamped configuration while the theory was based on simple support boundary 
conditions. The fairly good agreement between the experimental and theoret- 
ical results lends support to the implicit assumption that the modal density 
of sandwich panels, like single- layered panels, is relatively independent of 
the boundary conditions. 


CONCLUSIONS 


Theoretical estimates were obtained for the modal densities of sandwich 
beams and flat or cylindrically curved sandwich panels. The relative impor- 
tance of transverse shear flexibility and orthotropic shear moduli of the core, 
bending stiffness of the faces, rotary inertia, and panel curvature as they 
affect modal density was evaluated. Experimental values of modal density were 
obtained from resonance tests of flat rectangular panels having orthotropic 
cores. From the results of the investigation the following conclusions are 
noted : 

1. Failure to account for the transverse shear flexibility of the core 
can lead to a significant underestimation of the modal density. 

2. The effect of rotary inertia is generally negligible compared to the 
effect of transverse shear flexibility. 
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3. For many practical sandwich configurations and frequency ranges, the 
effect of the face bending stiffness can be neglected. 

4. For flat rectangular panels having orthotropic cores, a 90° rotation 
of the core material with respect to the faces has no effect on the average 
modal density. 

5. Where face bending effects are small, the results obtained for flat 
and cylindrically curved sandwich panels having isotropic cores can be applied 
to panels having moderately orthotropic cores, such as honeycomb, by introduc- 
ing an effective shear modulus and an effective curvature parameter. The 
effective shear modulus is simply the geometric mean of the two-face parallel 
shear moduli. 

6. The agreement between modal densities predicted by theory and modal 
densities determined from experiment was generally good except at frequencies 
near that of the fundamental mode where the theory is not applicable. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Dec. 2, 1969 
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APPENDIX A 


PANEL STIFFNESS AND INERTIAL PROPERTIES 



I 


M 

* d l> d 2 , and e are given in symbols list and were obtained from equation (18) of reference 22. 
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APPENDIX B 

DERIVATION OF EQUATION (1) 


The differential equations governing the panel vibrations are obtained 
from the small deflection theory of reference 11 by adding transverse and 
rotary inertia terms. 


3Q^ 9Q V 
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3y 
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3x 4 


M 


3fw 
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. 3 2 Q 

1 + y x 
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y 
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= 0 


= o y (bi) 


= o 




In equations (Bl) , x and y specify the coordinates of a point in the 
panel’s middle surface (fig. 1) and t denotes time. The operator V -4 is 
defined by V _4 (V 4 w) = V 4 (V" 4 w) = w where V 4 = 3 4 /3x 4 + 23 4 /3x 2 3y 2 +3 4 /3y 4 . 

The motion described by these equations is that a straight line perpen- 
dicular to the undeformed middle surface (z = 0) remains straight and of con- 
stant length after deformation but not necessarily perpendicular to the 
deformed middle surface. This inclination in the x (or y) direction from a 
right angle is the average shear angle Qx/Dq x (or Qy/Dq ) produced by the 

resultant transverse shear force Q x (or Qy) per unit width. 

For simply supported edges parallel to the x axis at which the support 
is applied over the entire thickness, the boundary conditions are (ref. 23) 


W = My = p- 


= 0 


(B2) 
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where the moment 


My 


(acting about the 


x axis) is given by 




i 3( V + / L. 9Q x \ 

\ 37 \3X 2 " D Q X 9X / 


(B3) 


The boundary conditions along the edges parallel to the y axis are obtained 
by interchanging y and x. 

Expressions for the lateral deflection and shear angles that satisfy the 
boundary conditions are 


, ... nnrx . nny iwt 

w(x,y,t) = A' sm sin ~~ e 

3. D 


x f . B f mTTX . niry icot 

D— Dn - C0S ~~a~ Sin b 6 


y r ... C' . mux mry ioit 
— (x,y,t) = ^ — sm -j- cos e 


y 


y 


r ( B4 ) 


where m and n are integers designating the number of sinusoidal halfwaves 
in the x and y directions, respectively, and w is the panel frequency 
(rad/ sec) . 


The differential equations (Bl) are also satisfied by the above forms for 
w, Q x /Dq , and Qy/Dq provided that 


(#’(*)* - (Z ) 2 
(*)[»*- »(£)’] 

" < (/) ] 


a/b 


-{1 + r v 


(*)’• -«(£)’ 


- ( i T J1 )(i7b)"x 


1 + r,. 


‘ ( 1 T JL )(i7b) nr y 


= o 


(B5) 


Expanding the determinant yields equation (1) . 
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APPENDIX C 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International Systems of Units (SI) was adopted by the Eleventh 
General Conference on Weights and Measures, Paris, October 1960, in Resolution 
No. 12 (ref. 24). Conversion factors for the units used herein are given in 
the following table: 


Physical quantity 

U.S. 

customary 

unit 

Conversion 

factor 

(*) 

SI Unit 

Length 

in. 

0.0254 

Meters (m) 

Force 

lbf 

4.448 

Newtons (N) 

Moment 

in. - lb f 

0. 113 

Meter-Newtons (m-N) 

Density 

lbf-sec 2 
in. 4 

1.069xl0 7 

Kilograms per cubic meter (kg/m 3 ) 

Modulus 

lbf 

in. 2 

6.895xl0 3 

Newtons per square meter (N/m 2 ) 


*Multiply value given in U.S. Customary Unit by conversion factor to obtain 
equivalent value in SI unit. 


Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

tera 

(T) 

10 12 

giga 

(G) 

10 9 

mega 

(M) 

10 6 

kilo 

00 

10 3 

centi 

(C) 

10 -2 

mil li 

(m) 

10" 3 
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Figure 




m 

a/b 


Figure 2.- Modes bounded by a constant frequency curve 





Figure 3.- Effect of shear flexibility and rotary inertia on the modal 
density of a flat sandwich panel; isotropic core; face bending 
stiffness neglected. 



Figure 4.- Effect of shear flexibility and face bending stiffness on 
the modal density of a flat sandwich panel; isotropic core; rotary 
inertia neglected. 
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Figure 5.- Effect of shear flexibility and rotary inertia on the modal 
density of a sandwich beam; face bending stiffness neglected. 



Figure 6.- Effect of shear flexibility and face bending stiffness on 
the modal density of a sandwich beam; rotary inertia neglected. 
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c 


Figure 7.- Modal densities of cylindrical ly curved sandwich panels; 
isotropic core; face bending stiffness neglected. 



Figure 8.- Asymptotic behavior of cylindrically curved sandwich panel 
modal density; isotropic core; face bending stiffness and rotary 
inertia neglected. 
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(a) rC = 3.0 
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(b) rC = 5.0 


Figure 9.- Combined effect of curvature and face bending stiffness on sandwich 
panel modal density; isotropic core; rotary inertia neglected; tf/h c = 0.10. 
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(c) rC = 10.0 
Figure 9.- Concluded. 



Figure 10.- Effect of core orthotropy on cumulative number of modes 
for a cylindrically curved panel as predicted by equation (26); 
rotary inertia and face bending stiffness neglected; a/b = 1.27. 


1 000 1 — ■ 



(a) Exact; equation (26). 



(b) Empirical; equation (27). 

Figure 11.- Comparison of cumulative number of modes for sandwich 
panels as predicted by exact and empirical frequency equations; 
rotary inertia and face bending stiffness neglected; a/b = 2.55, 
y = 0.5. 






Figure 13.- Combined effect of orthotropic core shear moduli and face 
bending stiffness on the modal density of a flat sandwich panel; 
rotary inertia neglected. 



PANEL I 



PANEL 2 



Figure 14.- Core geometries and orientations. 

exaggerated. ) 


(Core cell size is 



A-41046 


(b) Mode (4, 4) - 1780 Hz. 

Figure 15.- Some Chladni figures on panel 4 showing modes (m, n) 



( 



A-41050 

(c) Mode (5, 6) - 2630 Hz. 



A-41051 

(d) Mode (8, 6) - 3240 Hz. 


Figure 15.- Concluded. 
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EQUATION (8) 


EQUATION (8) 



(a) Panel 1. 



(c) Panel 3. 

Figure 16.- Comparison of theoretical 



(b) Panel 2. 



(d) Panel 4. 

and experimental modal densities. 
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